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摘要 :【 目 的] 意大利 蝗 Calliptamus italicus 是 新 疆 荒漠 、 半 荒漠 草原 优势 危害 种 类 之 一 ,其 生殖 相 
关 基 因 信 息 缺 乏 。 本 研究 对 意大利 蝗 精 业 和 卵 业 组 织 进 行 高 通 量 转录 组 测序 ,站 在 揭示 意大利 
蝗 转 录 组 特征 ,并 获得 与 性 腺 发 育 相 关 的 基因 数据 。【 方 法】 利用 Ilumina 高 通 量 测序 平台 
( HiSeg/MiSeq) 对 意大利 蝗 进行 了 转录 组 测序 ,并 进行 序列 分 析 . 【结果 】 获 得 718 872 条 
unigenes ,平均 长 度 631 bp, N50 为 1 186 bp。 通 过 同 源 性 比 对 ,成 功 注释 254 597 条 unigenes, 其 
中 ,Nr 数据 库 注 释 占 28.87% ,比例 最 高 。Nr 数据 库 注 释 的 unigenes 与 黑 褐 草 蚁 Lasius niger 的 相 
似 基因 序列 最 多 ,为 10.80% 。 与 GO 数据 库 比 对 发 现 ,意大利 凰 转录 组 unigene 可 分 为 3 大 类 : 涉 
及 分 子 功 能 的 unigenes 有 31 392 条 ,涉及 细胞 组 分 的 unigenes 有 20 586 条 ,与 生物 学 过 程 有 关 的 
mcs 有 57 0144, KEGG pathways 2-4 4&9] ,23 666 条 unigenes 归属 于 251 个 通路 。 涉 及 生殖 
系统 发 TERCER IA \ 胰 岛 素 信号 通路 、Wnt 信号 通路 、 促 性 腺 激素 释放 激素 信号 
通路 、 转 化 生长 因子 B 信号 通路 以 及 泛 素 介 时 的 蛋白 ee pb 卵 母 细胞 成 熟 . 昆虫 激素 
生物 合成 等 。 进 一 步 得 选 得 到 部 分 与 性 腺 发 育 相 关 的 基因 ,包括 ve, ver, sox 21b, testis-specific 
serine/threonine-protein kinase, spermatogenesis-associated. protein 和 vasa- Pw gene 等 。 【结论 】 本 研究 
获得 了 意大利 蝗 转 录 组 数据 及 性 腺 发 育 相 关 基 因 , 为 进一步 研究 意大利 蝗 生殖 调控 机 理 商定 了 分 
子 基础 。 
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Abstract: [^im] Calliptamus italicus is a major plague pest in the arid and semi-arid grasslands of 
Xinjiang, with a lack of information of its reproduction related genes. The objective of this study is to 


reveal the characteristics of the transcriptome of C. italicus and to obtain the gonada development related 
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genes in C. italicus through the transcriptome sequencing of its testicular and ovarian tissues. [ Methods] 
The transcriptome of C. italicus was sequenced using the Ilumina HiSeq/MiSeq platform and 
bioinformatically analyzed. [Results] A total of 718 872 unigenes were obtained, with a mean length of 
631 bp and an N50 of 1 186 bp. Through a similarity search, 254 597 unigenes were annotated. Most 
unigenes (28. 8796 ) were annotated to Nr database, and the unigenes of C. italicus had the highest 
similarity (10. 8096 ) with those of Lasius niger annotated in Nr database. According to GO database, all 
unigenes were broadly annotated into three categories; 31 392 unigenes to be related to molecular 
function, 20 586 unigenes related to cellular component, and 57 014 unigenes related to biological 
processes. In KEGG database, a total of 23 666 unigenes were assigned to 251 metabolic pathways, of 
which many pathways are involved in the reproductive system development, including oocyte meiosis, 
insulin signaling pathway, Wnt signaling pathway, GnRH signaling pathway, TGF- signaling pathway, 
ubiquitin mediated proteolysis, progesterone-mediated oocyte maturation, and insect hormone 
biosynthesis. By further screening and identification, gonadal development related genes of C. italicus, 
including genes of vg, vgr, sox 21b, testis-specific serine/threonine-protein kinase, spermatogenesis- 
associated protein and vasa-like gene, were obtained. [Conclusion] This study acquired the transcriptome 
data and gonadal development related genes of C. italicus, providing a molecular foundation for further 
studying the mechanism of reproductive regulation of C. italicus. 


Key words: Calliptamus italicus; high-throughput sequencing technology; gene annotation ; ovary; testis 














随 着 新 一 代 高 通 量 测序 技术 的 发 展 ,转录 组 分 
析 已 经 成 为 非 模式 昆虫 发 掘 功能 基因 的 有 将 手段 之 
一 ,广泛 应 用 于 遗传 进化 、 害 虫 防治 .生态 适应 等 研 
究 中 (Oppenheim et al., 2015; Zhang et al., 2016) 。 
近 些 年 ,高通 量 测序 技术 被 大 量 应 用 于 蝗虫 转录 组 
学 研究 中 (Km et al., 2016; Qiu et al., 2016, 
2017) 。 金 永 玲 等 (2015 ) 对 大 垫 尖 翅 蝗 Epacromius 
coerulipes 进行 了 转录 组 分 析 , 共 得 到 了 63 033 条 
unigenes ,筛选 得 到 与 杀 虫 剂 抗 药性 相关 的 功能 基 
及 潜在 SSR 位 点 。 张 硕 (2015 ) 通过 构建 亚洲 小 车 
I Oedaleus asiaticus 上 典 雄 成 虫 触角 转录 组 文库 , 鉴 
定 出 编码 气味 分 子 结合 蛋白 基因 ,并 对 其 表达 谱 进 
行 测定 分 析 ,为 进一步 研究 亚洲 小 车 蝗 嗅觉 感受 的 
分 子 机 制 提供 数据 基础 。 

蝗虫 具有 较 强 的 生殖 适应 能 力 ,种 类 多 、 数 量 
大 分 布 广 ,是 危害 农 牧 业 生产 安全 的 世界 性 重大 害 
虫 ,其 生殖 调控 机 制 一 直 是 国内 外 研究 热点 ( 赵 旨 
MAHR, 2005) 。 蝗 虫 始 终 保持 较 高 的 产生 配子 
的 能 力 ,是 其 暴发 成 灾 的 重要 原因 ( 赵 半 等, 2008), 
性 腺 是 蝗虫 繁殖 的 重要 场所 , 精 业 和 卵巢 发 育成 熟 
并 行使 正常 的 生理 功能 ,才能 产生 有 效 精子 和 卵子 ， 
保障 蝗虫 种 群 繁 衔 ( 赵 盾 等 , 2008) 。 因 此 ,对 蝗虫 
性 腺 发 育 的 分 子 调控 机 制 开展 研究 ,有 利于 深入 了 
解 蝗虫 的 生殖 调控 机 理 , 为 害虫 的 综合 防治 提供 理 
论 支 持 及 数据 基础 。 转 录 组 技术 为 从 分 子 水 平 研究 
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昆虫 性 腺 发 育 和 配子 发 生 的 调控 机 制 提供 了 新 思 
路 。 杨 婧 (2013 ) £i 9j fa Wi Actractemorpha. sinensis 
雌雄 成 虫 转录 组 进行 分 析 , 筛 选 出 1 879 个 肉 雄 差 
异 基因 ,发 现 组 织 蛋 白 酶 在 肉 忠 体内 表达 上 调 , 而 
CYP3A 在 雄 虫 体内 表达 上 调 , 推 测 二 者 分 别 参与 了 
卵子 发 生 过 程 和 性 激素 生物 合成 过 程 ,为 深入 研究 
蝗虫 性 别 分 化 及 生殖 调控 机 理 提供 了 数据 支持 。 

ACRI Calliptamus italicus Ti A SH dé É A 
漠 草 原 的 优势 危害 种 类 , 随 着 全 球 气候 变化 ,意大利 
蝗灾 害 在 新 疆 频 繁 发 生 ,给 新 疆 的 畜牧 业经 济 、 生 态 
和 社会 造成 严重 的 损失 ( 薛 智 平等 , 2010) 。 目 前 ， 
对 意大利 蝗 的 研究 主要 集中 在 生态 分 布 . 生 长 发 育 、 
行为 特征 .生理 生化 、 分 子 生 物 学 及 防治 等 方面 ,而 
关于 意大利 蝗 生 殖 调控 方面 的 研究 相对 较 少 ( 王 瞪 
等 , 2010; 起 上 忠 伟 等 ,2013; Hunter et al., 2016; 
Popova et al., 2016; Ren et al., 2016; 王冬梅 等 ， 
2016; 张 晓 红 等 , 2016) 。 前 期 研究 发 现 , 高 温 对 意 
大 利 蝗 肉 虫 繁殖 力 有 显著 影响 , 随 温度 升 高 ,意大利 
蝗 卵 巢 发 育 和 卵子 发 生 进程 加 快 (向 敏 等 , 2017 ) 。 
然而 ,其 分 子 机 制 还 有 待 于 进一步 研究 。 意 大 利 蝗 
为 非 模式 生物 ,可 参考 的 相关 基因 、 和 蛋白 序列 较 少 ， 
因此 ,本 研究 借助 高 通 量 测序 技术 对 意大利 蝗 精 迷 
和 卵巢 转录 组 进行 分 析 , 以 期 获得 卵巢 发 育 相 关 的 
重要 功能 基因 ,为 从 分 子 水 平 上 阐明 意大利 蝗 -的 生 
殖 适 应 机 制 及 暴发 成 灾 规律 提供 科学 依据 。 
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1 材料 与 方法 


1.1 ftii 

意大利 蝗 于 2016 4E 7 月 采 自 昌吉 回族 自治 州 
333r. (43*54'N, 86?21'E, YFIR 1 310 m) , 带 回 
实验 室 后 置 于 人 工 气 候 箱 (浙江 托 普 GTOP-260B) 
内 ,在 27 x 1*C 、 相 对 湿度 45% +5% 、 光 周期 14L: 
10D 条 件 下 饲养 ,挑选 性 成 熟 期 (未 交配 ) 、 健 康 的 
雌雄 成 虫 ,分 别 取 卵 梨 和 精 全 后 迅速 投入 液 氮 冷 冻 ， 
并 于 -80% 保存 竺 用。 每 3 个 卵巢 或 精 集 混合 为 一 
个 检测 样品 , 卵 和 偶 和 精 介 各 共 取 3 个 检测 样品 。 
1.2 总 RNA 的 提取 与 检测 

采用 Trizol Reagent( Invitrogen ) 提取 总 RNA ,分 
别 使 用 Nanodrop2000 (IMPLEN, CA, USA) , Qubit 
RNA Assay Kit, Agilent 2100 ( Agilent Technologies , 
CA, USA) fll 1. 096 琼脂 糖 凝 胶 电 泳 检 测 总 RNA 样 
品 的 纯度 、 浓 度 和 完整 性 。 
1.3 cDNA 文库 构建 和 测序 

RNA 样品 检测 合格 后 ,用 带 有 Oligo(dT) 的 磁 
珠 富 集 mRNA, 随后 加 入 fragmentation buffer 将 
mRNA 打 断 成 短片 段 ,以 mRNA 为 模板 构建 测序 文 
库 。 文 库 构 建 完 成 后 ,使 用 Qubit2. O, Agilent 2100 









































1.5 基因 功能 注释 
使 用 BLAST 软件 将 unigene 序列 与 下 列 数据 库 
比 对 ,获得 基因 的 功能 注释 信息 。 数 据 库 分 别 为 : 
NCBI 蛋白 数据 库 (NCBI non-redundant protein 
sequences, Nr) , NCBI 核酸 序列 数据 库 (NCBI non- 
redundant nucleotide sequences, Nt), 4E H X ÑK 
(protein family, Pfam), EZ [E] US AE DS ( clusters of 
orthologous groups of proteins, KOG/COG ) , Swiss- 
Prot 数据 库 (a manually annotated and reviewed 
protein sequence database) , 京都 基因 与 基因 组 百科 
全 书 ( KEGG ortholog database, KO) , 基因 本 体 论 
( Gene Ontology, GO) 。 
1.6 性 腺 发 育 相关 基因 的 筛选 

根据 雌雄 样本 差异 表达 基因 的 功能 注释 信息 ， 
结合 已 有 相关 文献 报道 ,筛选 并 鉴定 与 性 腺 发 育 过 
程 相关 的 功能 基因 。 























2 结果 
2.1 意大利 蝗 转录 组 的 测序 和 组 装 








采用 Hlumina 测序 平台 对 意大利 蝗 转 录 组 进行 
测序 ,每 个 样本 生成 超过 12.7 Gb 的 高 质量 数据 
(NCBI 登录 号 : SRR6113309, SRR6113310), GC 





fll Q-PCR 方法 对 文库 进行 检测 ,以 保证 文库 质量 。 
最 后 使 用 Tllumina 高 通 量 测序 平台 (HiSeq/ MiSeq) 
进行 测序 。cDNA 文库 构建 与 测序 由 北京 康 普 森 生 
物 技 术 有 限 公 司 协助 完成 。 
1.4 转录 组 数据 的 拼接 与 组 装 

将 测序 得 到 的 原始 数据 (raw reads ) 经 过 处 理 可 
获得 高 质量 的 测序 数据 (clean data) 。 同 时 , 计算 
Q20 和 Q30, GC 含量 和 重复 序列 水 平 , 并 利用 
Trinity 软件 对 获得 高 质量 的 clean data 数据 进行 序 
列 拼接 .过滤 和 组 装 ,最终 获 得 高 质量 的 unigene。 








含量 在 41.92% ~47.25% 之 间 。 每 个 样本 Q20 均 
大 于 97% ,Q3 均 大 于 93% ( 表 1)。 测 序 结 果 准 确 
度 较 高 ,可 以 用 于 后 续 分 析 。 

利用 Trinity 软件 进行 序列 组 装 , 共 得 到 832 098 
个 转录 本 ,总 长 度 为 651 690 832 bp, N50 和 平均 长 
度 分 别 为 1 930 bp 和 783 bp。 对 转录 本 进行 聚 类 和 
组 装 分 析 , 得 到 718 872 条 unigenes，N50 和 平均 长 
度 分 别 为 1 186 bp 和 631 bp( 表 2) 。 转 录 组 长 度 分 
布 见 表 2, 其 中 长 度 小 于 600 bp 的 unigenes 所 占 比 
例 最 大 ,为 80. 481% 。 

















表 1 意大利 蝗 6 个 cDNA 样品 测序 数据 评估 统计 


Table 1 Summary of sequencing data of six cDNA samples of Calliptamus italicus 


数据 大 小 





样品 编号 。 原始 Reads 数 高 质量 Reads 数 CCAA) 020 Qo 
Size of clean bases 
Sample no. Number of raw reads Number of clean data (Gb) GC content (96) (96) 
;b 
Mos 52 459 346 52 270 701 15.68 41.25 97.49 93.58 
42 766 547 42 598 957 12.78 41.92 97.28 93.34 
Female locust 

43 370 123 43 219 365 12.96 42.20 97.25 93.26 

n 53 694 335 53 500 693 16.05 45.09 97.33 93.32 

1 52 711 023 52 505 418 15.75 44.33 97.27 93.21 

Male locust 

55 590 859 55 378 554 16.61 44.37 97.30 93.28 





Q20, Q30: 4j dés fie (H2 20 和 30 的 碱 基 所 占 百 分 比 Percentages of bases with the quality value =20 and 30, respectively. 
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R2 意大利 蝗 转 录 组 组 装 结果 


Table 2 Illumina transcriptome assembly results for Calliptamus italicus transcriptome 











长 度 范围 转录 本 Transcripts Unigenes 
Length range( bp) 数量 Number 百分比 Percentage 数量 Number 百分比 Percentage 
< -600 610 470 73.365 578 553 80.481 
600 —1 500 161 714 19.434 115 378 16.050 
1 500 -3 000 40 111 4.821 177 66 2.471 
3 000 -6 000 16 351 1.965 6 023 0.838 
6 000 - 15 000 3 388 0.407 1 131 0.157 
»15 000 64 0.008 21 0.003 
总 数 Total number 832 098 718 872 
总 长 度 Total length ( bp ) 651 690 832 453 313 494 
N50 长 度 N50 length( bp) 1 930 1 186 
平均 长 度 Mean length( bp) 783 631 


N50; Z xt 5096 所 有 核 苷 酸 的 最 大 Unigene KJE Length of the smallest contig which when added to a set of larger contigs yields at least 50% of the 


genome. 


2.2 基因 功能 注释 X3 不 同 数 据 库 中 意大利 蝗 转 录 组 unigenes 注释 
使 用 BLAST 软件 将 获得 的 unigenes 序列 与 Nr， 结果 统计 


Table3 Annonation result statistics of unigenes of 


























Nt, Pfam, KOG/COG, Swiss-Prot, KO 和 GO 各 数据 Calliptamus italicus transcriptome in different 

库 进 行 比 对 ,获得 意大利 星 转录 组 unigene 的 注释 信 databases 

息 。 在 组 装 获得 的 718 872 条 unigenes 中 ,成 功 注释 注释 数据 库 被 注释 的 基因 数量 注释 百分比 

254 597 条 unigenes, 占 总 unigenes 数 的 33. 4296 。 其 E E n o did rd 

中 ,Nr 数据 库 注 释 207 555 条 unigenes , j 28. 8796 , Nr 207 555 28.87 

比例 最 高 ,其 后 依次 是 KOG (15. 0196 ) , Swiss-Prot KOG 107 886 15.01 

(11. 1096) , Nt(7. 69%), GO (5. 9096), KEGG Swiss-Prot TUE HR 

(3. 29% ) 和 Pfam(0. 0396 ) (X 3) 2 Pn Mi 
9 d'un í i a GO 42 401 5.90 

Unigene 注释 到 Nr 数据 库 中 的 物种 分 布 见 图 KEGG 23 666 3.29 
1, 其 中 与 黑 褐 草 蚁 Lasius niger 的 相似 基因 序列 最 Pfam 189 0.03 
多 ,所 占 比例 为 10.80% , 其 后 依次 是 内 华 达 十 白蚁 ”于 0 2599507 Se 


物种 分 布 Species distribution 


2.1896 


2.1096 


i: 
25996 3,3596 — 4 sg, 





5.6196 





8.6496 


10.8096 
53.3696 















































B EIE GAB Cerapachys biroi B 麦 双 尾 是 Diuraphis noxia 
B jt Amyelois transitella 温带 臭虫 Cimex lectularius 
R ECUmi HA] P BLVolenhovia emeryi m HE Diaphorina citri 
c ri RE Acyrthosiphon pisum B 赤 拟 谷 盗 7ribolium castaneum 
国内 华 达 十 白蚁 Zootermopsis nevadensis B 黑 褐 草 蚁 Lasius niger 
四 其 他 Others 
图 1 意大利 蝗 unigene 在 Nr 数据 库 中 的 物种 分 布 图 





Fig. 1 Species distribution of unigenes of Calliptamus italicus in Nr database 
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Zootermopsis nevadensis (8. 64% ) 、 赤 拟 谷 盗 Tribolium 
castaneum ( 5. 6996 ) , Wü 豆 WP Acyrthosiphon pisum 
(5.61% ) HIRR BU Diaphorina citri(3. 5996) RR 
遍 胸 切 叶 蚁 Vollenhovia emeryi (3. 3596 ) ,其 他 相似 
性 大 于 2% Im) RR ur di Cimex lectularius 
(2.59% ) , op RE Amyelois transitella (2. 18% )、 麦 
双 尾 蚜 Diuraphis noxia (2. 10% ) , E& EC JH. ff ii hy. 
Cerapachys biroi(2. 09% ) ,其 他 物种 占 53. 3696 。 
2.3 Unigene 的 功能 分 类 

GO 是 国际 标准 化 的 基因 功能 分 类 系统 ,能 够 
将 基因 划分 为 3 大 类 ,分 别 是 分 子 功能 (molecular 
function) 、 细 胞 组 分 (cellular component ) 和 生物 过 程 
(biological process) 。 对 意大利 蝗 转 录 组 unigene 进 
fT GO 分 析 发 现 ,涉及 分 子 功能 的 unigenes 有 31 392 
条 ,其 中 分 子 功 能 (molecular function ) 、 催 化 活性 
( catalytic activity ) 和 结合 (binding ) 的 unigenes 数量 
最 多 ,分 别 是 4 455, 2 762 和 2 553 条 , 其余 均 在 
1 500 条 以 下 ;涉及 细胞 组 分 的 unigenes 有 20 586 
条 ,其 中 细胞 组 分 (cellular component)、 细 胞 膜 
(membrane) 、 细 胞 (cell) 和 细胞 部 分 (cell part ) 的 
unigenes 数量 最 多 ,分 别 是 2 909, 1 657, 1 502 和 
1 502 条 ,其 余 大 部 分 在 1500 条 以 下 ;与 生物 过 程 有 
关 的 unigenes 有 57 014 条 ,其 中 生物 过 程 (biological 
process) 代谢 过 程 (metabolic process) 、 细 胞 过 程 
(cellular. process )、 有 机 物 代 谢 过 程 (organic 
substance metabolic process ) 的 unigenes 均 在 2 000 
条 以 上 ,分 别 为 3 629, 2 483, 2 381 和 2 022 条 ,其 
余 均 在 2 000 条 以 下 (图 2)。 

此 外 , 所 有 Unigenes 经 过 KOG 数据 库 功能 预测 
和 分 类 ,共有 94 262 条 基因 被 分 为 26 类 (图 3)。 其 
中 只 有 一 般 功 能 (general function prediction only ) 的 
unigenes 数量 最 多 ,为 46 367 条 (49. 19% ) ,是 最 大 
类 群 ; 其 次 分 别 是 参与 信号 转 导 机 制 (signal 
transduction mechanisms ) 的 5 430 条 (5.76% ) ,参与 
细胞 周期 调控 细胞 分 裂 和 染色 体 分 区 (cell cycle 
control, cell division, chromosome partitioning ) 的 
5 4247& (5.7596 ) ,参与 翻译 后 修饰 、 蛋 白质 转换 、 分 
F fE fH (posttranslational 
turnover, chaperones) 的 4 600 条 (4.88% ) , & Ej je 
色 质 结构 和 动力 学 (chromatin structure and 
dynamics) 的 3 909 条 (4. 1596), 参与 转录 
( transcription) 的 2 569 条 (2.73% ) ,参与 未 知 功能 
(function unknown) 的 2 383 条 (2.53% ) ,参与 糖 类 
运输 与 代谢 (carbohydrate transport and metabolism ) 





















































modification, protein 


的 2 275 2& (2.4196) ,参与 翻译 ,核糖 体 结构 和 生物 
合成 (translation, ribosomal structure and biogenesis) 
的 2 202 条 (2. 34%), 参与 复制 .重组 和 修复 
(replication, recombination and repair) 的 2 151 条 
(2. 2896), 参与 细胞 内 运输 .分 泌 和 膜 泡 运输 
( intracellular trafficking, 
transport) 的 2 060 条 (2. 1996) , unigenes 数量 在 
1 000 ~2 000 之 间 的 有 氨基 酸 运 输 与 代谢 (amino 
acid transport and metabolism ) 、 能 量 的 产生 和 转换 
(energy production and conversion ) , RNA 加 工 和 修 
4f (RNA processing and modification )、 脂 类 运输 与 代 
谢 (lipid transport and metabolism) 次 级 代谢 产物 的 
生物 合成 、 转 运 和 分 解 (secondary metabolites 
biosynthesis, transport and catabolism )、 细 胞 骨架 
(cytoskeleton) 无 机 离子 转运 和 代谢 (inorganic ion 
transport and metabolism ) 和 细胞 壁 / 膜 / 包 腊 的 生物 
RÆ (cell wall/membrane/envelope biogenesis) ,分 别 
为 1917, 1838, 1808, 1764, 1 711, 1583, 1131 
和 1 049 条 。 其 他 的 unigenes 数量 均 少 于 1 000 条 ， 
其 中 与 核 结 构 (nuclear structure ) 和 细胞 动力 (cell 
motility) 有 关 的 unigenes 更 少 , 所 占 比 例 分 别 为 
0.21% 710.07% 。 

使 用 KEGG 注释 系统 进行 unigene 代谢 途径 分 
析 ,结果 表 明 ,23 666 条 unigenes 归属 于 251 条 通 
路 ,其 中 碳 代谢 (carbon metabolism) 通路 最 多 , 为 
468 条 ;含有 300 ~ 400 条 unigenes 的 通路 有 :PI3K- 
Akt 信号 途径 (PI3K-Akt signaling pathway ) ,367 条 ，; 
氨基 酸 合成 (biosynthesis of amino acids) ,349 条 ; 含 
有 200 ~ 300 条 unigenes 的 通路 有 : 喀 啶 代谢 
( pyrimidine metabolism ) ,299 Z& ; YE 93 4H RE WE f il 
(starch and sucrose metabolism ) , 295 条 ; 粘着 斑 
(focal adhesion ) , 281 条 ; 7 Wi 4% ( phagosome ) 253 
条 ;细胞 周期 (cell cycle) , 244. 条 ; 糖 酵 解 途径 
( glycolysis /gluconeogenesis ) , 243 条 ; RNA 运转 
(RNA transport) , 237 条 ; 丙酮 酸 代 谢 (pyruvate 
metabolism), 232 4&; AMPK 信和 号 途径 (AMPK 
signaling pathway ) ,226 条 ;内 质 网 中 的 蛋白 质 合成 
(protein processing in endoplasmic reticulum ) , 224 
Z& ; cAMP 信号 途径 (cAMP signaling pathway ) , 219 
条 ; cGMP-PKG 信号 途径 (cGMP-PKG signaling 
pathway), 218 条 ; 卵 母 细胞 减 数 分 裂 ( oocyte 
meiosis) ,216 条 ; MAPK 信号 通道 (MAPK signaling 
pathway) ,208 条 。 其 他 通路 的 unigenes 数量 均 在 
200 条 以 下 ,其 中 涉及 性 腺 发 育 的 通路 有 : 胰岛 素 信 





secretion and vesicular 
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分 子 功能 Molecular function 

催化 活性 Catalytic activity 

结合 Binding 

有 机 环 状 化 合 物 结合 Organic cyclic compound binding 
杂 环 化 合 物 结合 Heterocyclic compound binding 
转移 酶 活性 Transferase activity 

水 解 酶 活性 Hydrolase activity 

小 分 子 结合 Small molecule binding 

离子 结合 Ion binding 

阳离子 结合 Cation binding 

碳水 化 合 物 衍 生物 结合 Carbohydrate derivative binding 
金属 离子 结合 Metal ion binding 

HTA Nucleoside binding 

EIZH A Purine nucleoside binding 

EZIZ iA Purine ribonucleoside binding 
细胞 组 分 Cellular component 








分 子 功能 Molecular function 








细胞 膜 Membrane 

细胞 Cell 

E 细胞 部 分 Cell part 
胞 内 Intracellular 

E 胞 内 部 分 Intracellular part 

9 细胞 膜 部 分 Membrane part 

z 细胞 膜 内 组 分 mtrinsic component of membrane 
i 细胞 器 Organelle 
胞 内 细胞 器 Intracellular organelle 
ER 膜 结合 细胞 器 Membrane-bounded organelle 
g 胞 内 膜 结合 细胞 器 Intracellular membrane-bounded organelle 
E 大 分 子 复 合 物 Macromolecular complex 


细胞 器 部 分 Organelle part 
胞 内 细胞 器 部 分 Intracellular organelle part 
生物 过 程 Biological process 
代谢 过 程 Metabolic process 

细胞 过 程 Cellular process 

有 机 物 代谢 过 程 Organic substance metabolic process 
主要 代谢 过 程 Primary metabolic process 
细胞 代谢 过 程 Cellular metabolic process 
大 分 子 代谢 过 程 Macromolecule metabolic process 
单 生物 过 程 Single-organism process 
细胞 大 分 子 代 谢 过 程 Cellular macromolecule metabolic process 
单一 生物 细胞 过 程 Single-organism cellular process 
含 氮 化 合 物 代谢 过 程 Nitrogen compound metabolic process 

细胞 含 氮 化 合 物 代谢 过 程 Cellular nitrogen compound metabolic process 
细胞 芳香 族 化 合 物 代谢 过 程 Cellular aromatic compound metabolic process 
有 机 环 状 化 合 物 代谢 过 程 Organic cyclic compound metabolic process 

杂 环 化 合 物 代谢 过 程 Heterocycle metabolic process 

含 碱 基 化 合 物 代谢 Nucleobase-containing compound metabolic 


蛋白 质 代谢 过 程 Protein metabolic process 
核酸 代谢 过 程 Nucleic acid metabolic process 

定位 Localization 

单一 生物 代谢 过 程 Single-organism metabolic process 


建立 定位 Establishment of localization 




















生物 过 程 Biological process 











基因 表达 Gene expression 
生物 合成 过 程 Biosynthetic process 
细胞 内 蛋白 质 代谢 过 程 Cellular protein metabolic process 
运输 Transport 




















图 2 GO 


1500 3000 4500 
GO 注释 数量 Number of GO annotations 


Fig. 2 Gene Ontology ( GO) categories 


号 通路 (insulin signaling pathway) ,185 条 ; Wnt 信和 号 
通路 (Wnt signaling pathway) ,183 条 ; 促 性 腺 激素 释 
放 激 素 信 号 通路 ( GnRH signaling pathway ) ,115 条 ; 
转化 生长 B 信号 通路 (TCF-B signaling 
pathway) ,101 条 ; 泛 素 介 导 的 蛋白 水 解 (ubiquitin 
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mediated proteolysis) ,76 条 ;黄体 酮 介 导 的 卵 母 细胞 


成 就 ( progesterone-mediated oocyte maturation ) , 66 


条 ; 昆虫 激素 生物 合成 (insect hormone 
biosynthesis) , 5 条 ( 表 4) 。 
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(A)RNA processing and modification 

B)Chromatin structure and dynamics 

C)Energy production and conversion 

D)Cell cycle control, cell division, chromosome partitioning 
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S)Function unknown 

T)Signal transduction mechanisms 
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Fig. 3  Karyotic Ortholog Groups ( KOG) classification 


X4 


意大利 蝗 转录 组 序列 KEGG 通路 统计 


Table4 Statistics of KEGG pathways for Calliptamus italicus transcriptome 


代谢 通路 


Pathways 








碳 代 谢 Carbon metabolism 

PI3K-Akt 信号 途径 PIBK-Akt signaling pathway 
氨基 酸 合成 Biosynthesis of amino acids 

WEE Pyrimidine metabolism 

淀粉 和 芒 糖 代谢 Starch and sucrose metabolism 

粘着 斑 Focal adhesion 

RIR Phagosome 

细胞 周期 Cell cycle 

糖 酵 解 途径 Glycolysis/Gluconeogenesis 

RNA 运转 RNA transport 

丙酮 酸 代谢 Pyruvate metabolism 

AMPK 信和 号 途径 AMPK signaling pathway 

j 质 网 中 的 蛋白 质 合成 Protein processing in endoplasmic reticulum 
cAMP 信和 号 途径 cAMP signaling pathway 

cGMP-PKG 信号 途径 cGMP-PKG signaling pathway 
STAR CLA Oocyte meiosis 

MAPK 信号 通道 MAPK signaling pathway 

其 他 通路 Others 

胰岛 素 信和 号 通路 Insulin signaling pathway 

Wnt 信号 通路 Wnt signaling pathway 

促 性 腺 激素 释放 激素 信号 通路 Gn RH. signaling pathway 
转化 生长 因子 B 信号 通路 TGF-B signaling pathway 
泛 素 介 导 的 蛋白 水 解 Ubiquitin mediated proteolysis 
黄体 酮 介 导 的 卵 母 细胞 成 熟 Progesterone-mediated oocyte maturation 
昆虫 激素 生物 合成 Insect hormone biosynthesis 

通路 总 数 Total 






































通路 ID 基因 数量 
Pathway ID Number of genes 

ko01200 468 
ko04151 367 
ko01230 349 
ko00240 299 
ko00500 295 
ko04510 281 
ko04145 253 
ko04110 244 
ko00010 243 
ko03013 237 
ko00620 232 
ko04152 226 
ko04141 224 
ko04024 219 
ko04022 218 
ko04114 216 
ko04010 208 

= <200 
ko04910 185 
ko04310 183 
ko04912 115 
ko04350 101 
ko04120 76 
ko04914 66 
ko00981 5 

- 23 666 
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2.4 意大利 蝗 性 腺 发 育 相关 基因 的 注释 上 调 , 且 10 824 个 基因 仅 在 雌 虫 表达 。 基 于 差异 表 


对 意大利 蝗 雌 雄 成 虫 基因 表达 量 进行 分 析 ， 
筛选 得 到 24 653 个 差异 表达 基因 ,其 中 12 954 个 
虫 表达 上 调 ,11 699 个 基因 在 


基因 在 雄 








达 基 因 的 功能 注释 , 结 
x2 ds rm ia» 
HORS). 








fE rh RIX 








结合 已 有 的 相关 文献 报道 ， 
及 性 腺 发 育 过 程 的 重要 功能 



















































































m5 意大利 蝗 转录 组 中 与 性 腺 发 育 相关 的 基因 
Table 5 Gonadal development related genes in Calliptamus italicus transcriptome 
基因 物种 长 度 (bp) E f 一 致 性 (% ) 
Genes Matched species Length E value Identity 
性 腺 发 育 Gonadal development 
vg ACH Lethocerus deyrollei 4 967 0 31.199 
vgr 美洲 大 归 Periplaneta americana 5 451 0 45.163 
sox 21b 黑 脉 金 斑 蝶 Danaus plexippus 350 6.27E -24 97.778 
PREDICTED;  testis-specific | serine/threonine-protein 
松 叶 蜂 Neodiprion lecontei 4 035 0 84.848 
kinase 2 
PREDICTED; testis-specific serine/threonine-protein 
松 叶 蜂 Neodiprion lecontei 1 501 2.08E -150 68.667 
kinase 4-like 
Spermatogenesis-associated protein 5 内 华 达 古 白蚁 Zootermopsis nevadensis 2 703 0 80.658 
Spermatogenesis-associated protein 22 内 华 达 古 白蚁 Zootermopsis nevadensis 2 289 3.77E -29 44.526 
Vasa-like gene 沙漠 蝗 Schistocerca gregaria 3 792 0 90. 167 
昆虫 激素 生物 合成 Insect hormone biosynthesis 
Dehydrogenase/reductase SDR family member 11 内 华 达 古 白蚁 Zootermopsis nevadensis 874 2.36E — 68 52.792 
Juvenile hormone acid methyltransferase 沙漠 蝗 Schistocerca gregaria 2 210 3.22E -54 50.296 
Juvenile hormone epoxide hydrolase 1 内 华 达 古 白蚁 Zootermopsis nevadensis 1 025 1.26E - 154 69.072 
Juvenile hormone esterase 黄 星 桑 天 牛 Psacothea hilaris 2 406 9.84E -88 33.789 
PREDICTED : cytochrome P450 307a7 赤 拟 谷 盗 Tribolium castaneum 4 345 2.93E - 109 67.557 
黄体 酮 介 导 的 卵 母 细胞 成 熟 Progesterone-mediated oocyte maturation 
Membrane-associated tyrosine- and  threonine-specific 
内 华 达 古 白蚁 Zootermopsis nevadensis 2 402 1.86E - 77 48.084 
cdc2-inhibitory kinase 
Mitotic | checkpoint — serine/threonine-protein kinase 
内 华 达 古 白蚁 Zootermopsis nevadensis 3 883 6.96E -77 45.912 
BUBI 
Phosphatidylinositol-4 , ^^ 5-bisphosphate 3-kinase 
内 华 达 古 白蚁 Zootermopsis nevadensis 4 606 0 80. 1678 
catalytic subunit delta isoform 
PREDICTED: cytoplasmic polyadenylation element- 
TRINE Nasonia vitripennis 2 073 1.14E -09 100 
binding protein 2 isoform X4 
PREDICTED: LOW QUALITY PROTEIN.: adenylate 
大 蜜蜂 Apis dorsata 5 733 0 57.837 
cyclase type 9-like 
促 性 腺 激素 释放 激素 信号 通路 GnRH signaling pathway 
Adenylate cyclase type 2 短 管 赤 有 眼 蜂 Trichogramma pretiosum 4372 3.09E - 174 80.958 
Calmodulin 柑橘 凤 蝶 Papilio xuthus 2 603 2.97E -95 98.675 
Early growth response 1 东亚 飞 蝗 Locusta migratoria 3 990 0 98.288 
Epidermal growth factor receptor 内 华 达 十 白蚁 Zootermopsis nevadensis 7 265 0 66.482 
Mitogen-activated protein kinase kinase kinase 1 内 华 达 古 白蚁 Zootermopsis nevadensis 4 628 0 41.745 
PREDICTED: gonadotropin-releasing hormone 工 T 
JEE IT RE Athalia rosae 1 016 8.13E -21 75.385 
receptor-like 
PREDICTED: LOW QUALITY PROTEIN: adenylate 
大 蜜蜂 Apis dorsata 5 733 0 57.837 
cyclase type 9-like 
Transcription factor AP-1 内 华 达 古 白蚁 Zootermopsis nevadensis 2 104 1.13E -35 69.792 
Tyrosine-protein kinase Srcó4 B 内 华 达 十 白蚁 Zootermopsis nevadensis 2 558 0 81.887 
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续 表 5 Table 5 continued 

基因 物种 长 度 (bp) E fH 一 致 性 (%% ) 

Genes Matched species Length E value Identity 
卵 母 细 胞 减 数 分 裂 Oocyte meiosis 
Aurora kinase A 袋 儿 Sarcophilus harrisii 4 575 1.85E -27 72. 824 
Calmodulin 柑橘 凤 蝶 Papilio xuthus 2 603 2.97E - 95 98.675 
Cell division cycle protein 20-like protein 内 华 达 十 白蚁 Zootermopsis nevadensis 2 765 5.69E -71 34.459 
G1 / S-specific cyclin-E 内 华 达 十 白蚁 Zootermopsis nevadensis 3 106 1.18E -132 53.102 
Membrane-associated tyrosine- and  threonine-specific 
ER 内 华 达 十 白蚁 Zootermopsis nevadensis 2 402 1.86E -77 48.084 
Mitotic | checkpoint — serine/threonine-protein kinase 
Em 内 华 达 十 白蚁 Zootermopsis nevadensis 3 883 6.96E -77 45.912 
PREDICTED; cytoplasmic polyadenylation element- E 
Modo 丽 晶 晴 集 金 小 蜂 Nasonia vitripennis 2 073 1.14E -09 100 
PREDICTED; LOW QUALITY PROTEIN; adenylate 
oom 大 蜜蜂 Apis dorsata 5 733 0 57.831 
Protein phosphatase 2A regulatory subunit A , putative 点 蜂 缘 晴 Riptortus pedestris 282 1.65E -40 75.269 
Serine/threonine-protein phosphatase PP1 内 华 达 十 白蚁 Zootermopsis nevadensis 4 443 4.67E-178 89.007 
细胞 周期 Cell cycle 
CDKS regulatory subunit-associated protein 1 JE IE Athalia rosae 2.356 0 66. 160 
Cell division cycle 2-like protein kinase 内 华 达 十 白蚁 Zootermopsis nevadensis 1735 0 60. 405 
Cell division cycle 5-like protein JKXH B ABUS] Cephus cinctus 5244 0 76.903 
Cyclin dependent kinase 1 红火 蚊 Solenopsis invicta 6 048 2.48E - 17 82.645 
Cyclin dependent kinase 2 PERN Drosophila. suzukii 2 388 6.64E -08 77.143 
MAPK 信号 通路 MAPK signal pathway 
14-3-3 zeta IARI Gesonula punctifrons 4 122 0 83.967 
14-3-3 protein epsilon Vb ida Schistocerca gregaria 6 158 3.54E -167 94.465 
Mitogen-activated protein. kinase 1-like PÆ Limulus polyphemus 2 043 4.31E -120 74.901 
Mitogen-activated protein kinase kinase kinase 15 EN BEBE Harpegnathos saltator 5 571 0 75.481 
泛 素 介 导 的 蛋白 水 解 Ubiquitin mediated proteolysis 
E3 ubiquitin-protein ligase SIAHI PERENG Lepisosteus oculatus 6 831 1.36E-169 81.744 
Polyubiquitin-C-like 橄榄 果 蝇 Bactrocera oleae 260 3.40E - 13 88.571 
Ubiquitin carboxyl-terminal. hydrolase isozyme L5 美洲 东部 能 蜂 Bombus impatiens 4 700 2.51E-11 71.551 


3 讨论 


对 于 缺少 基因 组 信息 的 非 模式 生物 ,通过 转录 
组 测序 ,可 以 从 特定 组 织 中 获得 参与 特定 生物 学 过 
程 的 功能 基因 信息 (Kiinstner et al., 2010; Jung et 
al., 2011; Teaniniuraitemoana et al., 2014) 。 本 研究 
选取 性 成 熟 期 (交配 前 期 ) X ACR IG PE ZH ER VETT 
转录 组 测序 和 分 析 ,以 期 获得 意大利 蝗 生 殖 相关 基 
因 信息 , 共 获 得 718 872 条 高 质量 unigenes, 其 中 有 





254 597 条 unigenes 在 Nr, KOG, KEGG, Swiss-Prot 
和 GO 等 数据 库 中 被 成 功 注释 。KEGG 分 析 表 明 ， 
23 666 条 unigenes 归属 于 251 条 代谢 通路 ,包括 一 
些 在 性 腺 发 育 和 成 熟 过 程 中 必要 的 信和 号 通路 ( Meng 
et al., 2015) ,如 与 内 分 泌 系 统 (endocrine system ) 相 
关 的 胰岛 素 信号 通路 、 促 性 腺 激素 释放 激素 信和 号 通 
路 黄体 酮 介 导 的 卵 母 细胞 成 熟 途径 , 与 信号 转 导 
(signal transduction ) 相关 的 PIBK-Akt 信号 途径 、 
AMPK 信号 途径 、cAMP 信号 途径 、cGMP-PKG 信号 
途径 .MAPK 信号 通道 、Wnt 信号 通路 、 转 化 生长 因 
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子 B 信号 通路 等 。 同 时 , 亦 鉴定 到 参与 信号 转 导 的 
相关 基因 ,如 腺 苷 酸 环 化 酶 2( adenylate cyclase type 
2) 8535] (calmodulin) 磷脂 酰 肌 醇 4,5- 二 磷酸 激 
Fg 4E, NV. d ( phosphatidylinositol-4 , 5-bisphosphate 3- 
kinase catalytic subunit delta isoform) 、 栈 氨 酸 蛋白 激 
Hi src64b ( tyrosine-protein kinase Src64B ) 基因 等 ,这 
些 基 因 在 意大利 蝗 性 腺 组 织 中 表现 出 较 高 的 表达 
量 。 相 关 代 谢 通路 和 基因 的 鉴定 表明 , 内分泌 调节 
系统 和 信号 转 导 系统 在 意大利 蝗 性 腺 发 育 和 功能 调 
节 中 具有 重要 作用 。 通 过 生物 信息 学 分 析 , 有 助 于 从 
大 规模 转录 组 数据 库 中 识别 潜在 的 功能 基因 ,为 意 大 
利 蝗 生 和 殖 调 控 机 制 研究 提供 丰富 的 基因 信息 资源 。 
本 研究 筛选 得 到 24 653 个 雌雄 差异 表达 基因 , 
这 些 基 因为 研究 意大利 蝗 生 殖 调 控 机 制 提 供 了 重要 
WREE. MEEA MITRE ERE ME k y 
力 的 关键 因素 ,是 种 群 繁衍 的 重要 基础 。 卵 黄 发 生 
是 卵子 发 生 中 的 关键 时 期 ,是 昆虫 生殖 调控 的 核心 ， 
卵 母 细胞 中 卵黄 蛋白 的 迅速 积累 , 可 促进 卵 母 细胞 
和 卵巢 的 发 育成 熟 , 亦 可 为 胚胎 发 育 提 供 营养 物质 ， 
有 利于 提高 孵化 率 及 幼虫 生长 发 育 , 对 昆虫 种 群 数 
量 有 显著 影响 (Muhammad et al., 2014) 。 卵 黄 发 生 
受 多 种 因素 的 调控 ,其 中 卵黄 蛋白 (vitellogenin, Vg) 
及 其 受 体 (vitellogenin receptor, VgR ) 基因 的 表达 是 
卵黄 发 生 的 基础 ,对 昆虫 卵巢 成 熟 起 着 至 关 重 要 的 
作用 ,也 是 研究 控制 害虫 的 潜在 训 标 ( Muhammad et 
al., 2014; Roy-Zokan et al., 2015; Upadhyay et al., 
2016)。 本 研究 中 ,Vg 和 VgR 在 意大利 星 卵巢 转录 
组 中 表达 上 调 , 揭示 其 在 卵 梨 发 育 中 的 重要 作用 。 
在 意大利 蝗 转录 组 中 亦 发 现 大 量 与 精子 发 生 相关 的 
基因 ,如 sox 2715、 举 九 特 异性 丝氨酸 / 苏 氨 酸 蛋白 激 
Hi ( testis-specific serine/threonine-protein kinases ) 、 精 
子 发 生 相 关 和 蛋白 C spermatogenesis-associated. protein ) 
基因 等 ,这 些 基 因 在 精 集 组 织 中 高 表达 ,提示 其 参与 
了 精 虽 分 化 和 精子 发 生 的 调控 。vasa 基因 是 
DEAD-box 家 族 基因 的 一 员 , 广 泛 存 在 于 所 有 真 核 
生物 和 大 多 数 原核 生物 中 ,是 生殖 系统 发 育 的 重要 
调控 因子 (Cordin et al., 2006), vasa 基因 突变 将 导 
致 生殖 细胞 分 化 和 卵细胞 发 育 的 异常 ,影响 正常 的 
卵子 发 生 进 程 , 造 成 雌性 个 体 不 育 (Styhler et al., 
1998) 。 研 究 表 明 , vasa [X TE Vb W WR Schistocerca 
gregaria Ji SIBI ER P RES EA vasa mRNA 在 早期 
卵子 发 生 中 高 表达 , 而 在 发 育成 熟 的 卵 中 不 表达 
( Chang et al., 2002), Z& BESE TE zx ACRI E f 5 I 
SUE ZH rS AZ IE vasa 基因 , 且 表 达 量 较 高 ,揭示 





































































































vasa 基因 在 意大利 蝗 性 腺 发 育 中 发 挥 重 要 作用 。 此 
外 , 亦 获得 参与 细胞 周期 . 泛 素 介 导 的 蛋白 水 解 等 通 
路 的 相关 基因 ,研究 表明 ,这 些 信 和 号 通路 均 参 与 一 系 
列 生殖 调控 活动 ,包括 配子 发 生 、 性 腺 发 育 等 (高 
iB, 2013; Meng et al., 2015; Peng et al., 2015; 
Yang et al., 2016) 。 

蝗虫 性 腺 发 育 还 受到 激素 等 物质 的 调控 , 其 
中 , 保 幼 激素 (juvenile hormone, JH) 对 于 卵黄 发 生 
及 卵 母 细胞 成 熟 具 有 重要 的 调节 作用 ( Guo et al., 
2014 ) 。 本 研究 在 意大利 蝗 性 腺 转录 组 中 注释 得 
到 参与 保 幼 激素 合成 和 代谢 的 功能 基因 ,如 保 幼 
激素 酸 甲 基 转 移 酶 ( Juvenile hormone acid 
methyltransferase )、 保 幼 激素 环 氧化 物 水 解 酶 1 
(juvenile hormone epoxide hydrolase 1) 、 保 幼 激 素 酯 酶 
(juvenile hormone esterase ) 基 等 ,这 些 基因 的 鉴定 
和 功能 分 析 将 在 分 子 水 平 上 促进 意大利 蝗 生 殖 内 分 
泌 研 究 。 

本 研究 利用 多 种 数据 库 对 所 获得 的 unigene 进 
ÍT BLASTX 比 对 分 析 ,成功 注释 254 597 条 unigene, 
并 筛选 得 到 参与 意大利 蝗 性 腺 发 育 和 配子 发 生 的 相 
关 功 能 基因 ,这 些 基 因 在 意大利 蝗 生 和 殖 调控 中 的 确 
切 作 用 还 有 待 于 进一步 研究 。 然 而 ,由 于 unigenes 
整体 长 度 较 短 .昆虫 基因 信息 仍 相 对 缺乏 等 原因 , 仍 
有 64.58% (464 275 条 unigene) 的 unigene 没有 获 
得 注释 信息 ,其 中 可 能 包含 新 的 与 性 腺 发 育 和 生殖 
调控 相关 的 功能 基因 , 这 些 基 因 的 鉴定 和 功能 分 析 
还 有 待 于 进一步 研究 。 
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